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Abstract 

We study the problem of determining the distribution of vertices of a par- 
ticular given type in the set of all Feynman tree graphs in quantum field 
theories. We show that in almost all Gaussian distribution arises 

asymptotically, and we compute the mean and variance of this distribution 
for several theories. We show the distribution's 'fine structure', arising from 
topological sum rules, can be obtained. 
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1 Introduction 



Recently, computational prowess in particle phenomenology has progressed 
to the evaluation of multi-particle scattering amplitudes, at least at the tree 
level, with ten or more external legs. The concomitant enormous number 
of Feynman graphs (e.g. 10.5 million for gluonic 2^8 scattering) suggests 
the study of the Feynman graphs themselves as an ensemble of combinato- 
rial objects. The number of graphs for various processes in self-interacting 
theories has been studied in a number of publications |2| , and also results 
for theories involving different particles have been found 0. It is therefore 
natural to examine various statistical aspects of these ensembles more closely. 
In the present paper we deal with the problem of determining the frequecny 
of occurrence of particular vertices. For instance, one may wish to find the 
frequency of occurrence of four-gluon vertices in the ensemble of Feynman 
graphs for gluonic QCD: it is this kind of question that we address in this 
paper. The layout of the paper is as follows. First, we show that the fre- 
quency distribution of any kind of vertex tends to a Gaussian distribution as 
the number of external lines becomes asymptotically large. We then proceed 
to evaluate the parameters of these distributions for several theories, notably 
gluonic QCD and theories with all possible vertices present. These results 
hold for the 'coarse-grained' point of view in which the number of vertices is 
interpreted as a continuous, rather than a discrete, fraction of the number of 
legs. Subsequently, we discuss how the 'fine structure' of the distribution can 
be determined: this fine structure arises from the topological conservation 
laws that govern the number of vertices in any tree amplitude with a given 
number of legs. Finally, we turn to the more realistic case of QCD, in which 
both quark and gluon fields occur. 



2 Gaussian limits for vertex occurrence 

Consider a self-interacting theory of a single field tp. Since we will only be 
counting diagrams and vertices, the particular types of interactions are only 
important in terms of the number of legs involved in each vertex. We shall 
start with a quick review of the counting of diagrams. Let the theory have 
vertices of type (p q+1 for some set of numbers q > 2 (that is, we allow any 
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combination of 3,4,5,. . . vertices). We define the 'potential' 



W( 9 ) 



q>2 



(1) 



where e q = 1 if the vertex occurs in the action (or effective action), otherwise 
it is zero. If we denote the tree-level 1 — > n amplitude by a(n), we can 
define the amplitude-generating function by 



[x 



^2 ^ro(n) 



n>l 



III 



(2) 



Then, simple diagrammatic arguments 
equation for the number of diagrams is 



|T| show that the Schwinger-Dyson 



(f)(x) = x + W(</>(x)) 



(3) 



Iteration of this relation, where <fi{x) is a power series in x, then gives the 
number of graphs for any desired value of n; alternatively, the analytical 
structure of d>( function of x informs us about the asymptotic behaviour 

of a(n) with n. In particular, it is clear that in the general case <p(x) will be 
singular in x, and moreover that this singularity will be a branch cut rather 
than a pole. The condition for the singularity is, of course 



d<p(x) 



dx 



dx 



oo 







W'((f)) = 1 



(4) 



Since q > 2 there are in general several solutions 0o to this condition, and 
the radius of convergence of <fi{x) as a power series is given by that solution 
4>o which leads to the smallest absolute value of x = <p — W^o)- Taking this 
value of 4>o, we see that, asymptotically, the number of graphs goes as 



a{n) ~ 



(0o - W{cf> Q )y 



(5) 



multiplied by sub leading terms: in fact, the more precise form has been 
derived in fill to read 



IV. 



1 



a(n) ~ 



0o ~ W^q) 

^(^o-%))1 27riy"(0 o ) 



(6) 



throughout this paper, the 'amplitude' is defined as the number of diagrams 
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but for the moment the leading form is sufficient. 

We now turn to the counting of vertices. Let us focus on the vertex of 
type ip q+1 for some q with e q — 1. We can count the occurrence of this vertex 
by giving it a weight z in the Schwinger-Dyson equation (|3|), which then 
becomes 

z) = x + W{<f>{x, z)) + {z- i)^£l! (7) 

The amplitude-generating function now depends on z as well as on x. By 
iteration of this implicit equation we can find 

4>(x, z) = —a(n; z) (8) 

n>l ^' 

to desired order in x. In the coefficient 

a(n;z) = £ b(n;m)z m (9) 

m>0 

which is a finite polynomial in z, the number b(n; m) is then the number 
of 1 — > n graphs with precisely m vertices of type ip q+1 ; and a(n; 1) is the 
total number of graphs. It is more instructive to consider not the number 
of graphs themselves but rather the probability to pick a diagram with m 
(q + l)-vertices at random out of all tree graphs in the 1 — > n amplitude. If 
we denote this probability by ir(n; m), we have simply 

P n (.) = £7r(n;m)^ = ^ • (10) 

m>0 a \ n i L ) 

Putting in the leading asymptotic behaviour for a(n) of Eq.(||) we therefore 
find that 

5>(n;m)s"~ ( ^T^ Y ' Ffaz) = <p-W(<p)-(z-l)£ {U) 



m>0 



F{4> Q {z)-z) 



and 4>o(z) is of course the value of <p(x; z) at the singularity nearest to the 
origin in the complex x plane: 

^(z))-^- 1)^^ = 1 . (12) 
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The distribution of the m values is determined once we know its moment- 
generating (or characteristic) function. It is simply related to the result we 
have obtained: 



k>0 h - m,k>0 

■^"^■(W (13) 

But, by elementary arguments from probability theory, this implies that m 
is distributed as the sum of n independent, identically distributed random 
variables. The central limit theorem therefore asserts (except in special cases 
to which we shall come later) that for large n, the values of m are 
normally distributed: 

" (nvm) ~ vsbj exp (~^w G " " (,) ) 2 ) (14) 

with the parameters fx(q) and o~(q) 2 still to be determined (note that since 
m is always integer, it is rather the ratio m/n that follows the continuous 
normal distribution). 



3 Results for various theories 

Having established that the m values are normally distributed, we need only 
to compute the mean and variance. We have, for the first two factorial 
moments, 

(m) = P>(l) , (m(m - 1)> = P'^l) . (15) 
Some straightforward algebra then leads us to 

fq 

fi(q + l) = 



qKf-W(f)) ' 



a {q+ if = wmhnf)) ' (16) 

where / = 0o(l), hence W'(f) = 1. Note that, since all the Taylor coeffi- 
cients of W(ip) are nonnegative, the numbers f,f — W(f) and W"(f) are all 
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positive. As a further check, we recall the topological relation, valid for any 
given tree diagram: 

J2(k - 2)m k = n - 1 (17) 

fc>3 

where is the number of vertices of type <p k in that diagram. Indeed, it is 
trivially seen that 

£ e q (q-l)fji(q + l) = l . (18) 

q>2 

In addition, we can straightforwardly extend our discussion to the case where 
we consider the combined distribution of two different vertex types, with q± + l 
and g2 + l legs (qi ^ 92), by introducing two counting weights Z\^. The result 
for the expected value of the product m ?1+ im 92+ i is, then, 

n 2 _|_ n n 

(m qi+1 m q2+1 ) = (/ _ w{f))2 qi]q2] - W»(f)(f w (f))( qi -l)\(q 2 -l\) ' 

(19) 

which provides the additional check 

o = E e 9 (<7- 1 )M<? + i) 2 

+ E W?i - - 1) f Kl+img2+l) - M«i + + 1)) (20) 

We now turn to a few concrete cases. In the first place, there is gluonic 
QCD, with three- and four-point vertices, so that 

W(<p) = \p 2 + ^ 3 • (21) 

For z — 1, the nearest singularity is at x = \/3 — A/3, reached for / = 0o(l) = 
+ We find that 

M 4)=m^ . . ( 4)^h^. (22) 

In figures 1 and 2 we give the frequency distribution of the four-point ver- 
tices for this theory, for n = 25 (2.3 10 32 diagrams) and n = 55 (5.5 10 91 
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Figure 1: Frequency distribution of the number of 4-point vertices in <p 3 + y? 4 
theory, for n = 25. The solid line is the asymptotic estimate. 
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Figure 2: Frequency distribution of the number of 4-point vertices in (p 
theory, for n = 55. The solid line is the asymptotic estimate. 
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Figure 3: Average fraction (m)/n for four-point vertices, as a function of n, 
in ip 3 + ip 4 theory. 
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diagrams). In figure 3, we give the actual value of the average number of 
four-point vertices as a fraction of n, for 1 < n < 80. The asymptotic value, 
0.163986. . ., is reached from below (as already evident from figures 1 and 2): 
the correction term appears to go as l/4n in this range. 

The second case of interest is a theory (like an effective theory, after 
tadpole and mass renormalization) where all types of vertices occur, that is, 
e q — 1 for all q. We have 

W((p) =e v -tp-l , (23) 
so that / = log(2), / - W(f) = 21og(2) - 1, and W"(f) = 2, and we have 

(log 2)9 



Mg + i) 



q\ (2 log 2 - 1) 



1+fi- 



(log2) 



5-2 



2(9 - 1) 



(24) 



In figure 4 we plot log /i(q + 1) for the first few values of q for this theory. 
The values of a(q + l) 2 follow quite closely: the numerical values of cr(6) 2 
and /x(6) already differ by less than one percent. 



4 Fine structure 

In figure 5, we give the actual and asymptotic distribution of 3-point vertices 
in ip 3 + ip 4 theory, for n = 35. Every Feynman diagram contains an even 
number of 3-point vertices. This follows of course from the topological sum 
rule for this theory: 

m 3 + 2m 4 =n — l (25) 

so that m 3 is either always even or always odd, depending on the parity of 
n. The asymptotic estimate does not reflect this; we shall now describe how 
also this 'fine structure' can be obtained, even in the asymptotic limit, as a 
consequence of the algebraic structure of the action. In order to be slightly 
more general, let us assume that there are two types of vertices present, one 
with q + 1 legs and one with p + 1 legs, and that we count the number of 
vertices of the first type. The Schwinger- Dyson equation now reads 

dfl d)P 

<f) = x + z\ + 26 
q\ p\ 
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Figure 4: log /J,(q) as a function of q for a theory with all possible vertices. 



and the condition for the singularity reads 

2 fe^ + (^ = 1 - (27) 

After having obtained a(n; z), we can of course extract b(n; m) by standard 
means: 

6(n;m = / v \ ' dz , 28 

V ' ; 2m J z m+1 ' V ; 

where the integral is over an infinitesimal counterclockwise loop around z — 0. 
Now notice that, for z — 0, the singularity condition has not one but p — 1 
distinct solutions, which we shall denote by 

fc = (sO^exp^^) , fc = 0,l,2,...,s-l , s = p-l . (29) 

The corresponding values for x = 4>k — 4>k/P^ a ^ nave the same absolute value 
(at z — 0, to be sure) and contribute equally to the asymptotic result. It 
therefore makes more sense to perform the integral not in terms of z but 
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Figure 5: Frequency distribution of 3-point vertices in tp 3 + y? 4 theory, for 
n = 35. The solid line is the coarse-grained asymptotic estimate. 
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in terms of the 0&, using Eq.(]5|). To this end, we write, for nonzero z, the 
singularity condition as 



(g-1)! {, (f> s \ "H 
s 



i 1 — —r = 7 i (30) 



where by "if" (0 s ) we denote an unspecified function of s . Inserting this 
expression for z everywhere, we find 



The above form for 6(n; m) is now a sum of loop integrals: 

n! /■ 1 I -FQM dg 



(31) 



s-l 



fc=0 



n— m(q— 1) 

fe 



(2mk\ l ' n+m{q ' 1) , 



owing to the s-fold symmetry of the unspecified function "H" (0 s ). Now, the 
sum of the powers of the roots of unity in the above gives zero except when 
1 — n + m(q — 1) happens to be a multiple of s — p — 1, or in other words 

m(q — 1) + m'(p — 1) = n — 1 (33) 

for some m', which is again precisely the topological sum rule. In that case, 
we find exactly s times the result from a single singular point, so that the 
normalization of the asymptotic distribution is preserved. 

Two remarks are in order here. In the first place, the conclusion remains 
unchanged if other vertices are present, as long as these all have ks + 1 
legs (k any integer larger than 1). If this algebraic symmetry is destroyed 
by the presence of another vertex type, there will in general be only one 
dominant singular point, and the vertex frequency distribution will not show 
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any 'quantisation' anymore. In the second place, it is instructive to note that 
we arrive at the quantisation condition by working in the neighbourhood of 
z = 0: once we move out to z — 1, where the singular points are no longer 
symmetrically distributed, we move toward coarse-graining and a continuous 
distribution. 



5 Vertex counting in QCD 

We now turn to the more realistic case of QCD with both quarks and gluons. 
For simplicity we shall only consider a single quark flavour, as extensions 
to the case of more flavours is straightforward. We shall count the various 
vertices by assigning to the qqg vertex a weight z q , and to the three- and 
four-gluon vertices weights z 3 and z 4 , respectively. In any 1 — > n amplitude 
we shall assign a factor u to an outgoing quark, a factor v to an outgoing 
antiquark, and a factor x to an outgoing gluon as before. The amplitude- 
generating function for an incoming quark is denoted by ip, that for an in- 
coming antiquark by ip, and that for an incoming gluon by 0, again as before. 
The QCD Feynman rules now lead to the following coupled Schwinger-Dyson 
equations at the tree level: 

i> = u + z q i>4> , ip — v + z q ip(/> , 

= x + |0 2 + |0 3 + z q W , (34) 
We can rewrite this as an equation in terms of alone: 

= ^+|0 2 + ^0 3 + n , Z = *v . (35) 

2 6 (1 - z q (j)Y 

The occurrence of the combination £ reflects fermion number conservation. 
The generating function is 



x n g pr. 

n g >0n q >0m 3 >0m4>0m q >0 U 9- n Q : 



= 0(x,e,i-)= ]T E E E E c(n g ,n q ;m 3 ,m 4 ,m q ) — ^-z^zT 4 z^ 

_ _ , n y 



(36) 

where z = (zs, z±, z q ), and c(n g , n q ; m 3 , m 4 , m q ) denotes the number of Feyn- 
man tree graphs with one incoming gluon, n g outgoing gluons, n q outgoing 
qq pairs, and precisely m 3 three-gluon vertices, m 4 four-gluon vertices and 
m q quark-gluon vertices. Again, explicit iteration of this equation gives <fi as 
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a multinomial in all its arguments, from which the number of graphs can be 
read off easily. Due to the larger number of counting variables the arriving at 
asymptotic estimates is more cumbersome in this case, but not qualitatively 
different. We have 



x 







(37) 



and the large- n g behaviour is determined as before by requiring dx/d(f) to 
vanish, which occurs when 



l-z 3 <P- f 2 



(38) 



For this value of £ we then have 

X(<f>,z) = F(cf>,E(<f),z),z) 



r iV vr 2 Zq 1 



2 ^ 



(39) 



By the same arguments as before, we find that the leading asymptotic be- 
haviour is given by 



g(n g ,n q ;z) 



c(n g ,n q ]m 3 ,m4;,m q )z™ 3 zf 4 z r q 
1 



«*3,4 



XI 



where the saddle point / is determined by 



d 



— ln g log X(<t>,z) +n q log 2(0,. 



(40) 



(41) 



From Eq. fl40D it follows, by the probabilistic argument used before, that 
for n g ,n q — > oo the numbers 777.3, WI4 and m q are all normally distributed. 
The mean values and variances are again found by taking the appropriate 
derivatives at z— (1, 1, 1). The saddle-point value is, in that case, given byf] 

R(p, f) = 6(1 - p) - 12/ + 3(1 + p)f + (3 - p)f = , (42) 



There are other roots, but since these do not depend on p they cannot be relevant 
ones. 
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where 

n q = pN , N = n g + 2n g . (43) 
The solution is given by 



/ = —— \ -1 -p - 26sin 



P 



- arcsin ( ^- (p 3 - 9p 2 + 15p - 23 
3 Vtr v 



6 = (l3-2p + p 2 ) 1/2 . (44) 



The solution / interpolates smoothly from — 1 + ~ 0.732 at p = to 
-3/5 + 14/5 sin(arcsin(282/343)/3) ~ 0.2854 at p = 1/2. The results for the 
mean values are 

(m q ) = N 



(m 3 ) = N 



(l-/)(-24/ + 5/ 3 + 9/ 2 + 6) 
2/(18/ - \2pf - 6f + 7pf 3 - 9f 2 + 9pf 2 - 6 + 6p) 

(_24/ + 5/3 + 9/ 2 + 6)(-2 + /2 + 2/) 
/ 2 (16/ - 8p/ - 5/ 3 + 5p/ 3 - 7/ 2 + 5p/ 2 - 6 + 6p) 



(m4) " N (-24/ + 5/3 + 9/ 2 + 6)(-2 + / 2 + 2/) " (45) 

In figure 6 we give the results for the expectation values (rrij)/N, with j = 
q, 3, 4. As a final check, note that the topological sum rule for the vertices in 
this case reads 

m 3 + 2m 4 + m q = n g + 2n q — 2~N. (46) 

This is borne out by the numerical results, but also the analytic estimates 
(145T) show that the combination 



(777.3 + 2?774 + ™>q — N) 

is precisely proportional to R(f,p) and hence vanishes at the saddle point. 
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Figure 6: Expected number of vertices of given type as a fraction of N = 
n g + 2n q in QCD with a single quark flavour. Along the horizontal axis we 
plot p = n q /N. The rising curve is the result for the quark-gluon vertices, the 
falling curves are for the four-gluon vertex (lowest curve) and the three-gluon 
vertex (middle curve). 
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